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Recent progress in the calculation of multi-loop, multi-scale diagrams is reviewed. 
Expansion techniques combined with new developments in Computer algebra allow 
to evaluate the R ratio for massive quarks up to order and, partly, even af . 
Similar techniques can be applied to Higgs- or Z-boson decays and mixed QCD 
and clectroweak interactions. 



1 Introduction 

During the past years amazing progress has been made in the experimental 
tests of the Standard Model of particle physics. Its electroweak sector has 
been scrutinized at LEP, mainly in Z decays and W pair production, and at 
the TEVATRON, mainly through the precise measurements of the W boson 
and top quark masses. Perturbative QCD has been tested at electron positron 
colliders at lower energies as well as in high energy experiments at LEP, at the 
TEVATRON through proton-antiproton collisions, and through lepton-nucleon 
scattering in particular at HERA. 

The large statistics collected in these experiments in conjunction with the 
small systematic uncertainty allows to test the theoretical predictions with high 
precision, requiring the inclusion of quantum corrections at least in one-loop 
approximation, occasionally two or even three-loop calculations are required. 
This applies on one hand to perturbative QCD, for electroweak observables on 
the other hand either purely weak two-loop corrections should be included or 
so-called "mixed" QCD and electroweak effects up to two or even three loops. 
Considerable progress has been made on the theoretical side which matches 
these requirements. Two-loop amplitudes, in particular those contributing to 
two-point functions can often be evaluated in closed form or through straight- 
forward numerical integrals — even with arbitrary mass assignments. Three- 
loop amplitudes, however, are at present only accessible in a few selected cases. 
These are, most notably, two-point functions with massless internal propaga- 
tors only and vacuum diagrams with massless and massive propagators — 
however, of identical mass throughout. Other cases of interest are "on-shell 
amplitudes" required for mass renormalization and certain types of integrals 
relevant for fermion pair production at threshold. 

Despite this seemingly limited set of diagrams a large variety of problems 
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can be solved with the help of expansion techniques. Frequently one finds a hi- 
erarchy of mass and/or energy scales which suggests to perform a Taylor series 
of the integrand, with the resulting integrals being considerably simpler. How- 
ever, the integrals are not necessarily analytic in the expansion xi^rameter, and 
more refinedlechniques, denoted hard mass, large momentumBo or threshold 
expansion B'dO are required. These can be formulated in a diagrammatic man- 
ner and allow to reduce a given amplitude into an (infinite) sum of products 
of simpler "master" amplitudes which are known in analytical form. In prac- 
tice often the first few terms provide already a sufficiently accurate numerical 
answer. In this approach, the demands on computation grow rapidly. Not 
only the Dirac algebra and the evaluation of integrals for individual diagrams 
has to be performed by computer algebra, with the appearance of hundreds if 
not thousands of diagrams also the generation of the basic diagrams, the ap- 
plication of the hard mass or large momentum procedure, the transformation 
of diagrams into algebraic expressions and the oveijall book keeping has to be 
performed automatically. (For a recent review seeQ.) 

In the following talk several characteristic examples will be given, partly 
taken from QCD, partly from electroweak interactions. The next chapter 
will be concerned with the 0(0^) evaluation of the cross section for elec- 
tron positron annihilation into massive quarks. Information on the vacuum 
polarization function at small and large /{^m?) as obtained via expansion 
techniques is combined with the two dominant terms close to threshold. Sub- 
sequently a variable transformation is applied as suggested by the analyticity 
structure of n(g^) in the cut complex plane. The numerical results will be con- 
trasted with those deduced from the threshold expansion. Techniques, results 
and limitations of the large momentum procedure will be presented in chap- 
ter ^. This includes in particular the calculation of the vacuum polarization 
of order a^, expanded in rr? jq^ up to fairly high power. The strategy for an 
evaluation of the absorptive part in order a^^m? jq^ and Q^^(m? jq^Y ^'^^ recent 
results are given in chapter ^ Results of relevance for electroweak measure- 
ments are discussed in chapter ^. This includes the top quark contribution 
to the p parameter in three-loop approximation, mixed corrections to the Z 
decay rate and "theory driven" results for the running QED coupling at scale 
Mz- Chapter ^ contains a brief summary and the conclusions. 

2 Three-Loop Heavy-Quark Vacuum Polarization 

The measurement of the total cross section for electron positron annihilation 
into hadrons allows for a unique test of perturbative QCD. The decay rate 
T{Z hadrons) provides one of the most precise determinations of the strong 
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coupling constant ag. In the high energy Hmit the quark masses can often 
be neglected. In this approximation QCD corrections to i? = a{e~^e~ 
hadrons)/cr(e'*'e~ — > fi^fJ-') have been calculated QiJ up to order a^. For pre- 
cision measurements the dominant mass corrections must be included through 
an expansionin vr? j s. Terms up to order a^^m? / s (seeu) and a^m'^/s^ (seeE3) 
and recently Ej even up to a^m'^/s'^ are available at present, providing an 
acceptable approximation from the high energy region down to intermediate 
energy values. For a number of measurements, however, the information on 
the complete mass dependence is desirable. This includes charm and bottom 
meson production above the resonance region, say above 4.5 GeV and 12 GeV, 
respectively, and, of course, top quark production at a future electron positron 
collider. To order as this calculation, was performed by Kallen and Sabry in 
the context of QED a long time agot3. With measurements of ever increasing 
precision, predictions in order aj. are needed for a reliable comparison between 
theory and experiment. Furthermore, when one tries to apply the 0{a) result 
to QCD, with its running coupling constant, the choice of scale becomes im- 
portant. In fact, the distinction between the two intrinsically different scales, 
the relative momentum versus the center of mass energy, is crucial for a stable 
numerical prediction. This information can be obtained from a full calculation 
to order only. Such a calculation then allows to predict the cross section in 
the complete energy region where perturbative QCD can be applied — from 
close to threshold up to high energies. It is then only the region very close 
to threshold, where the fixed order result remains inadequate and Coulomb 
resummation becomes important. InEj results for the cross section were cal- 
culated in order a^. They were obtained from the vacuum polarization Ii{q^) 
which was calculated up to three loops. The imaginary part of the "fermionic 
contribution" — derived from diagrams with a massless imark loop inserted in 
the gluon propagator — had been calculated earlier iuEJ. In this latter case 
all integrals could be performed to the end and the result was expressed in 
terms of polylogarithms. InE3 the calculation was extended to the full set of 
diagrams relevant for QCD. Instead of trying to perform the integrals analyti- 
cally, information of Tl{q^) from the large behaviour, the expansion around 
= and from threshold was incorporated. 

2.1 Outline of the Calculation^ 

The different behaviour at threshold makes it necessary to decompose 11 ac- 
cording to its colour structure. It is convenient to write: 

n(g2) ^ n(0) (q2) + ^AlUc^n^^) {q^)+( ^f^) ' n(2) (g^) + . . . , (1) 

TT \ TT J 
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n(2) = C|n^') + CaCfU^na + CpTniU^^^ + CfTU^P . (2) 

The same notation is adopted to the physical observable R{s) which is related 
to n(g2) by 

R{s) = 12nlmU{q^ = s + ie). (3) 
The contributions from diagrams with m light or one massive internal 

(2) (2) 

fermion loop are denoted by CFTnilll ' and CpTIVp ^ respectively. The purely 
gluonic corrections are proportional to or CaCf where the former are the 
only contributions in an Abelian theory and the latter are characteristic for 
the non- Abelian aspects of QCD. 

All steps described below have been performed separately for the first three 
contributions to II'^-'. In fact, new information is only obtained for 11^"^ and 
since Imllp'' is already known analytically 0. The contribution from a 
four-particle cut with threshold at 4m is given in terms of a two dimensional 
integralliiO which can be solved easily numerically, so 11^'' will not be treated. 

Let us now discuss the behaviour oiJl{q^) in the three different kinematical 
regions and the approximation method. 

Analysis of the high behaviour: The high energy behaviour of 11 provides 
important constraints on the complete answer. In the limit of small vr? j q^ the 
constant term and the one proportional to niLLf' (modulated by powers of 
ln(/x^/g^) have been calculated a long time agollj. The results for terms up to 
order (m^/q^)^ are described in chapter ^, provide an important cross check, 
however, they are not used for the moment. 

Threshold behaviour: General arguments based on the influence of Coulomb 
exchange close to threshold, combined with the information on the perturbative 
QCD potential and the running of dictate the singularities and the structure 
of the leading cuts close to threshold, that is for small v = ^/l — Am? / s. The 

term is directly related to the QED result with internal photon lines only. 
The leading \/v singularity and the constant term of Ra can be predicted from 
the nonrelativistic Greens function for the Coulomb potential and the 0(as) 
calculation. The next-to-leading term is determined_by the combination of 
one-loop results again with the Coulomb singularitiesLD. One finds 

<-=3g-3,'....). ,4) 

The contributions ^ CaCf and ~ CfTtii can be treated in parallel. For 
these colour structures the perturbative QCD potentialEJ 

VQCD(<f^) = -4^C^^^, (5) 
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will become important. The leading CaCf and CpTni term in R is propor- 
tional to In V and is responsible for the evolution of the coupling constant close 
to threshold. Also the constant term can be predicted by the observation, that 
the leading term in order as is induced by the potential. The 0{as) result 

2, \ 2u TT y 

is employed to predict the logarithmic and constant CpCa and CpTni terms 
of 0{a^) by replacing by av{4:p^ — v^s) as given in Eq. (|^). This implies 
the following threshold behaviour: 



,(2) 



RY> = 3^(3 - + .2) Q In ^ - 1 + . . . ) . (8) 



This ansatz suggested inlil canJae verified for the CpTni term where the 
result is known in analytical form O and it has also been confirmed for the 
NA term where the leading terms for small v have been derived recently (see 
below). Extending the ansatz from the behaviour of the imaginary part close 
to the branching point into the complex plane allows to predict the leading 
terms of n(q^) ^ \nv and ~ In^ v. 

Behaviour at = 0: Important information is contained in the Taylor 
series of n(g^) around zero. The calculation of the first seven nontrivial terms 
is based on the evaluation of three-loop tadpplp integrals with the help of 
the algebraic program MAT AD written in FORMEj which performs the traces, 
calculates the derivatives with respect to the external momenta. It reduces 
the large number of different integrals to one master integral and a few simple 
ones through an |elahorate set of recursion relations based on the integration- 
by-parts methodc3l2l. The result can be written in the form: 

n>0 ^ ' 

where the first seven moments are listed in El. 

Conformal mapping and Fade approximation: The vacuum polarization 
function II'^' is analytic in the complex plane cut from = Am^ to +oo. The 
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Taylor series in is thus convergent in the domain \q^\ < 4m^ only. The 
conformal mapping which corresponds to the variable transformation (z = 
gV(4m2)) 



transforms the cut complex z plane into the interior of the unit circle. The 
special points z = 0, 1, — oo correspond to tj = 0, 1, —1, respectively (Fig. ^1). 




Figure 1: Transformation between the z and lu plane 



The upper (lower) part of the cut is mapped onto the upper (lower) perime- 
ter of the circle. The Taylor series in ui thus converges in the interior of the 
unit ciixjk^ To obtain predictions for 11 (g^) at the boundary it has been sug- 
gestedE3'E3 to use the Fade approximation which converges towards ll{q^) even 
on the perimeter. To improve the accuracy the singular threshold behaviour 
and the large q^ behaviour is incorporated into simple analytical functions 
which are removed from II^^^ before the Fade approximation is performed. 
The quality of this procedure can be tested by comparing the prediction with 
the known result for ImllJ . 

The logarithmic singularities at threshold and large q^ are removed by 
subtraction, the 1/v singularity, whiph is present for the Cp terms only, by 
multiplication with v as suggested incj. The imaginary part of the remainder 
which is actually approximated by the Fade method is thus smooth in the 
entire circle, numerically small and vanishes at uj — I and uj = —1. 



2.2 Results 

After performing the Fade approximation for the smooth remainder with lu as 
natural variable, the transformation ( [l0| ) is inverted and the full vacuum polar- 
ization function reconstructed by reintroducing the threshold and high energy 
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Figure 2: Complete results (full line) are compared to the threshold approximations and 
the high energy approximations including the rri^ j s (dash-dotted) and the / (dashed) 
terms {x = 2m/y/s). 




terms. This procedure provides real and imaginary parts of II^^^ . Subsequently 
only the absorptive part of n*^^) (multiplied by 127r) will be presented. 

In the following it will be useful to plot the results as functions of a; = 
2m/y/s and alternatively of v = ^1 — Aw? / s. The first choice is particularly 
useful for investigations of the high energy region, the second one for energies 
close to threshold. Characteristic values of x, v and -^/s for charm (mc = 
1.6 GeV), bottom (m;, — 4.7 GeV) and top [mt = 175 GeV) quarks are listed 
m Table for easy comparisons. Energy values where a perturbative treatment 
is evidently inapplicable are denoted by dashes. 

In Fig. ^ the complete results are shown for = w? with i?^^^ and b!'^\ 
displayed separately. The solid line represents the full correction. The thresh- 
old approximation is given by the dashed curve. In the high energy region 
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0.97 
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0.995 


0.980 


0.954 


0.917 


0.866 


0.800 


0.714 


0.600 


0.436 


0.312 


0.243 


0.199 


y/s^ 


32 


16 


10.7 


8.0 


6.4 


5.3 


4.6 


4.0 












94 


47 


31.3 


23.5 


18.8 


15.7 


13.4 


11.8 


10.4 










3500 


1750 


1167 


875 


700 


583 


500 


438 


389 


368 


361 


357 



Table 1: Conversion between x and v and values for y/~s in GcV for charm, bottom and top 
production. 



the corrections containing the m? / s terms and the quartic approximations are 
included. It should be stressed that the latter are not incorporated into the 
construction of i?'^-* but they are evidently very well reproduced by the method 
presented here. This will be investigated in more detail in chapter ^. There 
it will be demonstrated that the high energy expansion (including sufficiently 
many terms) and the Pade result agree remarkably well between a; = and 
X = 0.7 which covers most of the perturbative region for charmed and bottom 
quarks. 

The results which include high moments up to Cg, C^ or even Cg, are re- 
markably stable down to very small values of w. Different Pade approximations 
of the same degree and approximants with a reduced number of parameters give 
rise to practically identical predictions, which could hardly be distinguished in 
Fig. ||. Minor variations are observed close to threshold, after subtracting the 
singular and constant parts. The r|emainder 5R for up to ten different Pade 
approximants is shown in Fig. ra. Inli3 it was demonstrated that there is a per- 
feet agreement for _Rj , -Rjvjj seems to converge to the solid line ([4/4], [5/3] 
and [3/5]) when more moments from small are included. The dashed lines 
are from the [3/3], [4/2], [2/4] and [3/4], the dotted ones from lower order Pade 
approximants. The dash-dotted curve is the [4/3] Pade approximant and has 
a pole very close to a; = 1(1.07. . .). For the Abelian part a classification of 
the different results can be seen: the dashed lines are [4/2] and [2/4], the solid 
ones [3/2], [2/3], [5/3], [3/5] and [5/4] Pade approximants. 

Recently additional terms froixLthe high energy expansion have been in- 
jected in tha^Pade approximation EiJS. The results, shown in Fig. ^ confirm 
those from 113 and demonstrate the stability of the Pade approximation. It 
must be stressed that a safe estimate of the remaining uncertainty in Ra and 
Rna amounts to less than 0.02 for v above 0.1 and is around 0.05 in the region 
V « 0.03. This region is, however, entirely dominated by the singular and con- 
stant terms with values around 50 and higher. The perturbative predictions 
for R are therefore under excellent control. It goes without saying that the 
function Il{q^) constructed this way and evaluated e.g. in the Euclidean region 
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could be a useful tool for other investigations of interest, for example for sum 
rules involving massive quarks. 

Quite some effort has been invested in the analytic evaluation of Ra and 
Rna close to threshold. As stated above, the singular and constajatj-tjerms 
(Eqs. (||), (0) and (||)) were derived through general considerations Il3't2l. To 
evaluate the vlnv and the v terms, hossfever, elaborate analytical calculations 
of the two-loop form factor were usedEa-cl with the results 

3 ,7r4 2 / 35 2 , 4 , 39 ] 

3 / y ,179 , 8 , 151 13 

While the result for Ra was still obtained in the framework of "classical" 
QED calculations, Rna was calculated ijaifLg a convenient technique which 
formalized the expansion for small v (see HJj'EI) . Pade and small v results are 
compared in Fig. |^, again subtracting first the singular and constant pieces. 
The slopes for very small v predicted by the two approaches are again in nice 
agreement, giving further credibility to the Pade result. However, it is also 
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Figure 5: R(s) plotted against -/s. The scale fj? = s has been adopted. The dashed curves 
correspond to the values Mc = 1.8 GeV, Mi, = 5.0 GeV and as{M^) = 0.121, whereas for 
the solid curves Mc = 1.4 GeV, M,, = 4.4 GeV and as{M|) = 0.115 is used. The dotted lines 
show a recent compilation of the available experimental data. The central curves correspond 
to the mean vski«es, upper and lower curves indicate the combined statistical and systematical 
errors. (From|i£|.) 

clear from this comparison that the small v expansion alone cannot lead to a 
reliable prediction over a larger energy range. This is explicitely demonstrated 
for the case of Ri where the analytic result is available. Terms at least of order 
are needed for a stable prediction. 

A compilation of theoretical results can be found inll3 where the prediction 
for massive quark production is compared with the measurements over a wide 
energy region (Fig. ^). 

Up to this point only the vector current correlator has been discussed. 
However, the techniques described above have also been applied to other cur- 
rentsEH: the axial vector, relevant e.g. for top production through the virtual 
Z boson, as well as scalar and pseudoscalar currents describing for example the 
decay of Higgs bosons into massive quarks. Recently also the singlet piece of 
top quark production through the axial current was evaluated e3, completing 
thus the 0(0^.) prediction of massive quark production. 

3 Large Momentum Expansions 

An alternative route towards an efficient evaluation of the polarization function 
is based on the large momentum expansion. In this case the polarization 
function Ii{q^) is expanded in powers of m'^/q^, multiplied by logarithms of 
rri^/q^, with the power of the logarithms, however, limited by the number of 
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Figure 6: Comparison of the small v expansion (solid curves) of SRa, SRna with the 
full result (solid and dotted curves) after subtraction of the singular and constant pieces. 
The difference between dotted and solid curve indicate the remaining uncertainty of the 
semianalytical result. For SRi successive approximations of the small v expansion (dashed 
curves) and exact result (solid) are shown. 



loops under consideration. Technically the expansion is given by a series of 
products of massless propagator and massive tadpole integrals. 

In principle the full information on the analytic function is contained in this 
series. The structure of the integrals is simplified by moving from two- to one- 
scale integrals. However, an enormous proliferation of the number of diagrams 
and the amount of algebraic calculation is observed, requiring the development 
of programs which implement the diagrammatic expansion, and translate the 
resulting diagrams into input files for other programs which in turn arc suited 
for the algebraic evaluation Jif individual diagrams. One example for such a 
"superprogpjn" is GEFICDMo which uses QGRAFl3 for the generation of dia- 
grams, LMPEiJ or EXPEiJ for the diagrammatic, expansion thcciugh the hard mass 
or large momentum procedure, and MATADE3 and MINCERc^l for the evaluation 
of diagrams. Even nested expansions with a hierarchy of several scales are pos- 
sible in this framework. (A more detailed description of the status of algebraic 
programs can be found incl.) After performing the expansion in j up to a 
given power, one may directly take the absorptive part and thus predict n((7^) 
in the high energy region. The comparison with the Padc result (discussed in 
chapter H) shows excellent agreement in the region of x between zero and 0.7 
and thus down to fairly low energies (Fig. 

One may even suspect that, given sufficiently many terms of the absorp- 
tive part alone, an approximation of arbitrary precision can be achieved. This 
is indeed true for the two-loop result and the "double bubble" contribution 
with massless quarks in the internal loop, which are available in analytical 
form and thus can be used as "toy models" . However, in both examples, the 
only threshold is located at ^ = 2m and convergence down to this point is 
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X X 

(2) (2) 
Figure 7: The Abelian contribution and the non-Abelian piece RffJ^^ Wide dots: no 

mass terms; dashed hnes: including successively more mass terms {m?/s)" up to n = 5; 
solid hne: including mass terms up to tm^/s)®; narrow dots: semi-analytical result. The 
scale fj? = ni? has been adopted. (FromEj.) 

naturally expected. In contrast, the other diagrams have four-particle cuts 
at y/s = 4m, suggesting convergence only above this point. This has been 
confirmed by a detailed studya of the double bubble contribution with equal 
masses in both fermion loops, where inclusion of an increasing number of terms 
does not lead to an improvement beyond x « 0.5. Given both real and imagi- 
nary parts, this problem could be and has been circumvented as discussed in 
chapter ^. However, below we will be interested in the situation where the 
absorptive piece only is available. 

Nevertheless, the relative smallness of the four-particle contribution and 
the slow opening of the phase space reduce this effect and a fairly good approx- 
imation of Ra and Rna is achieved even for x = 0.7, thus covering most of the 
interesting energy range (cf. Table |^). The same technkme has also been used 
for the axial, the scalar and the pseudoscalar currentaH. (This motivates the 
step to a^m^ in chapter ^.) The agreement of the two-loop result described 
in chapter |^ with the threshold and the large momentum expansion in the re- 
spective ranges of validity demonstrate that the perturbative NLO prediction 
of R{s) for massive quarks is under full control. 

4 Quartic Mass Terms in NNLO form Operator Product Expansion 

Fig. suggests that the first three terms provide an excellent approximation at 
X — 0.5 and are quite acceptable even at x = 0.7. Using this line of reasoning, 
a possible route for a NNLO prediction (orderi-p^ ) of R{s), includingj-jquark 
mass effects, is at hand. The massless result B and the-jm^/s terms □ have 
been obtained a long time ago. The strategy used inll3 allowed to predict 
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the a^m'^/s^ terms in R{s) by evaluating two-loop (!) tadpoles and massless 
propagators only. Additional ingredients are the operator product expansion 
and the renormalization group equations, plus certain anomalous dimensions. 
Using this method and algebraic programs it is thus possible to obtain the 
afm^/s^ terms from three-loop tadpoles and massless propag atorsBEl. Let 
me briefly describe this method: In a first step the OPE is applied to the time 
ordered product of two currents 

J dxe''>-T{j^{x)j,{Q)) = {q^q,~g^,q^)I^A{q\fi^a,)l+ (13) 

- 2 ^1 1 

Only three of the six operators On with dimension four are gauge invariant 
and contribute to physical matrix elements: 

Oi = G2, O^^mqq, Oe = fh\fi^) , (14) 

the others are required for the proper construction to the coefficient functions 
C„. For the NNLO calculation, Ci, C2 and Cg are required up to 0{a'^), 0{a^) 
and 0(as) respectively. To obtain these coefficient functions, only massless 
propagator type integrals, at most of three loop, are needed. To calculate the 
vacuum matrix elements of Oi and O2, massive tadpole integrals - at most 
of three loop - are needed. In a last step one uses renormalization group 
invariance of the dimension four part of Eq. (|l3|). Employing the anomalous 
dimension matrix lD of 01,2.6 one finally obtains the coefficients of the terms 
a^m'* In" q^ / jJ?' with n = 1,2, 3,|-Qiily. these terms contribute to the absorptive 
part and one finally arrives atfflOEa 

1.40923 



+ 104.833 x + x^ (139.014 - 4.83333 Q 

+ 541.753 x + x"^ (3523.81 - 158.311 - 
with X = Tfi'^{s)/s, Ix = ln(m2(s)/s) and Uf ~ 5. 

5 Expansion techniques and electroweak interactions 

Electroweak observables are frequently affected by the interplay between strong 
and electroweak interactions. Important examples are the hadronic contribu- 



i?^(s) =3<^ l-6a;•' + — 
' TT 



1 + 12 a;- 22x2 



- 12.7671 



9.66667^2) 



(15) 
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tions to the running of the QED couphng from the Thonipson hmit to Mz, 
QCD effects on the p parameter and related quantities, and "mixed" vertex 
corrections affecting for example the Z decay rate. Quark mass effects and 
their perturbative treatment are important for aQEoiMz). A detailed diiscpi- 
sion of the last topic is beyond the scope of this paper and can be found inli3c3. 
Let us present some aspects of the two remaining items. 

Gauge boson self- energies, the mass of the top quark and QCD: The in- 
direct determination of the top quark through quantum corrections prior to 
its observation in hadronic collisions can be considered one of the triumphs of 
the Standard Model. The experimental precision of the key observables, the 
masses of the top quark and the W boson together with the weak mixing an- 
gle as determined by asymmetry measurements has increased during the past 
years and this process will continue in the foreseeable future. In order to con- 
trol the influence of the top quark at an adequate level the inclusion of QCD 
corrections in the top and bottom induced self energies is mandatory. The 
dominant terms are characterized by the p parameter which, in lowest order, 
is given byEj 

Ap = 3^/2%!#. (16) 
In view of the large difference between pole and running mass at scale m* 

mt{pole) = fht{mt) (^1 + ^ ^ (17) 

inclusion of two and even three-loop QCD corrections to Ap is mandatory. 

Analytic r^ults are available in two-loop approximation not only for the 
leading termE2li£i Ap but also for all self energies, with arbitrary top and 
bottom masses EJ. The resulting shift in the prediction for Mw for fixed 
Gf, ctQED, Mz and nit = 175 GeV amounts to 68 MeV, well comparable 
to the present precision and significantly larger than the anticipated accuracy 
of roughly 30 MeV. To arrive at a precise prediction for the central value and 
to control the theoretical uncertainty, three-Jnop QCD contributions to Ap as 
well as to Ar are required. The p parameterc3 can again be expressed through 
diagrams with vanishing external momentum (vacuum or tadpole diagrams), 
the remaining quantities involve_two point functions at non- vanishing external 
momentum and can be obtainedcj through an expansion in the small mass ratio 
M|/m(. The leading three-loop term corresponds to a shift of —10.9 MeV. 
The first three terms, amounting to —13.7 MeV, are adequate for a prediction 
with an accuracy better than 1 MeV (Table ^). Conversely, this combined shift 
is equivalent to a reduction of the effective top mass by about 1.6 GeV. 
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5Mw in MeV 




«i 






611.9 


-61.3 


-10.9 


const. 


136.6 


-6.0 


-2.6 


1/M2 


-9.0 


-1.0 


-0.2 



Table 2: Theiehange in Mw separated according to powers in as and hit in the on-shell 
scheme. (Fromcj.) 



To exploit the precision expected from a future linear collider which will 
pin down rat to better than 200 or perhaps even 100 MeV, the inclusion of a1 
terms is evidently mandatory. 

Similar considerationsil are valid for the effective weak mixing angle which 
can be deduced from the left right asymmetry or the forward backward asym- 
metry in a straightforward way. 

An important issue in this connection is the size of uncertainties, arising 
either from uncalculatcd higher orders, or from the "parametric" uncertainties 
in as and rut- Shifts of Sag = 0.003 and 5mt = 5 GeV lead to changes in Mw of 
—2.4 MeV and 35 MeV respectively. The uncertainty from uncalculated higher 
orders is completely negligible. A VF-mass determination with a precision in 
the 10 MeV range should therefore be accompanied by a top mass measurement 
with a precision around or better than 1 GeV. 

Mixed QCD and electroweak vertex corrections: As stated above, gauge 
boson self energies, in particular those induced by fermion loops, give rise to 
the dominant radiative corrections for electroweak precision observables. Nev- 
ertheless, for a complete treatment of 0{aas) the inclusion of irreducible ver- 
tex corrections is necessary. These have to be distinguished from the reducible 
ones which are easily incorporated, if the electroweak result, including one-loop 
weak terms, is simply multiplied by the QCD correction factor {l + ag/T: + - ■ ■). 

The physics and the techniques of calculation are markedly difierent for 
vertices leading to light (m, d, s and c) quark pairs on one handa and for 
decays into bb on the other hand e3, a consequence of the presence of top 
quarks with their enhanced contribution to the vertex proportional rn^/M^. 
The irreducible one-loop vertex diagrams are shown in Fig. || 

To obtain all one-particle irreducible vertex diagrams in two-loop approx- 
imation , these have to be dressed with gluons in all conceivable ways. The 
resulting amplitudes are first studied for arbitrary by considering expan- 
sions in the ratio xz — q'^/M'^ and xw = q'^/itt-w i'^^ ^^r some diagrams, in 
^/xz,w)- Even for the limiting values xz = 1 and xw = M'^/My^ the exact 
results are well approximated, once sufficiently many terms are included. Us- 
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(a) (b) (c) (d) 

<y O ^ 



(f) (g) 



-~(D^ Kl>- "^D^ 



Figure 8: Diagrams contributing to ST^ in 0{a). Thin lines correspond to bottom quarks, 
thick lines to top quarks, dotted lines to Goldstone bosons and inner wavy lines represent 
W bosons. 

ing as = 0.12, a ^ 1/129, sinOw = 0.223 Mz ^ 91.19 GeV, it is foundS that 
the net effect of the nonfactorizable corrections is 



p(2-loop EW/QCD) _ ^p(l-loop EW) 

TT 



71 

-1.13(4) X 10-4 GcV for Z ^ a_a , . 

-1.60(6) X 10-4 GeV for Z ^ dd ^ ' 

The total change in the partial width T{Z s- hadrons) is obtained by summing 
over 2 down- type and 2 up- type quarks: 

Ar(Z -> u,d,s,c) = -0.55(3) MeV (19) 

which translates into the change of the strong coupling constant determined 
at LEP 1 equal to 

AT(Z hadrons) 0.55 „ 

Aas = -TT ; = TT « 0.001 (20 

r{Z hadrons) 1741 ^ ' 

This shift is somewhat smaller but still of the same order of magnitude as 
the experimental accuracy and should to be taken into account in the final 
analysis of LEP 1 data. Electroweak parameters extracted from Z decays are 
not affected by this correction. 

The two-loop corrections to the Zbb vertex are dominated by terms quadratic 
in ru^ The QCD corrections to these have been evaluated^lready some time 
agollj. The sub-leading terms cx Inm^ were obtained incEl and are of com- 
parable magnitude, thus indicating relatively slow convergence of the series. 
The copaplete evaluation is thus mandatory and has been recently been per- 
formedeZI. In contrast to the previous case one is confronted with three different 
scales, the masses of top, the W and the Z boson. 
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Using the hard mass procedure for ^ M^^M^,, one may factor out 
the mt-dependence. However, for a part of the diagrams one stiU is left with 
two-scale and even three-scale integrals involving and and £wz.M^, 
where is the electroweak gauge parameter which has been kept incZl. Al- 
though they appear to be only one-loop integrals, their exact evaluation up 
to ©(e) is inconvenient. Instead, the results lII were obtained by applying 
the hard mass procedure to these kinds of diagrams once more, this time us- 
ing ^wM'^,M'^ ^ M^. This seemingly unrealistic choice of scales can be 
well justified: It is not possible for an expansion to distinguish the inequality 
> Af| from 4M^ > M| or (mf -|- Mwf > M|, the latter ones being 
perfectly alright. The only matter is to perform the expansion on the appro- 
priate side of all thresholds, and here one is concerned with thresholds at 2Mw 
and at rtit + Mw- Therefore, the choice ^ is to be understood purely 
in this technical sense. Graphically this continued expansion looks as follows: 

where only those terms are displayed which are relevant in the discussion above 
and all others contributing to the hard mass procedure are merged into the 
ellipse. The thick plain line is the top quark, the thick wavy one a Goldstone 
boson with mass squared £,wMw^ ^'^^ example. The thin plain lines are b- 
quarks, the inner thin wavy lines are VF-bosons, the outer ones Z-bosons. The 
spring-line is a gluon. The mass hierarchy is assumed to be ^ ^w^w ^ 
^ Mf . The freedom in choosing the magnitude of S^w provides a welcome 
check of the routines and the results. 

The outcome of this procedure is a nested series: The coefficients of the 
Mvy/mt-expansion are in turn series in Mz/Mw- Note that in contrast to 
the decay into u, d, s, c there is no threshold at Mw which makes an additional 
expansion in Mw /Mz unnecessary. 

In view of this calculation the procedure of successive application of the 
hard mass procedure resp. the large momentum procedure has been imple- 
mented in a Fortran 90 program named EXP Ell. Therefore, given an arbitrary 
hierarchy of mass scales, the computation of a three-loop two-point function 
can now be done fully lautomatically. Even more, the link to the Fc^ojiman dia- 
gram generator QGRAFE2lin a common environment called GEFICOMlj allows to 
obtain the result of a whole physical process without any human interference 
except for specification of the process and final renormalization. 

Finally, the result for the VF-induced corrections to the Z-decay rate 
5T^ [Z — > bh) is conveniently presented in the form of the renormalization 
scheme independent difference to the decay rate into dd. Inserting the on-shell 
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top mass nit = 175 GeV, the Z-mass Mz = 91.19 GeV and sin^ Ow = 0.223 
gives cZl 



5r^{Z bb) - ST^iZ dd) = ro— -<j - 0.50 + (0.71 - 0.48) 



sm t)w ^ 



(0.08 - 0.29) + (-0.01 - 0.07) + (-0.007 - 0.006) + — 

TT 



1.16+ (1.21 



0.49) + (0.30 - 0.65) + (0.02 - 0.21 + 0.01) + (-0.01 - 0.04 + 0.004) 



-^0 



1 



sin^ e 



w 



0.50-0.07 



as 



1.16 + 0.13 



(21) 



where the factor T^a/ (tt sin^ 6w) with = Alza/ (4 sin'^ 9\y cos^ 6w) has been 
pulled out for convenience. The numbers after the first equality sign corre- 
spond to successively increasing orders in l/m^, where the brackets collect the 
corresponding constant, logm^ and, if present, log^ mt-terms. The numbers 
after the second equality sign represent the leading mj-term and the sum of 
the sub-leading ones. The 0{a) and 0(Q;as)-results are displayed separatebz. 
Comparison of this expansion of the one-loop terms to the exact result E3 
shows agreement up to 0.01% which gives quite some confidence in the aag— 
contribution. One can see that although the m^-, Wj- and log mt-terms 
are of the same order of magnitude, the final result is surprisingly well repre- 
sented by the TOj -term, since the sub-leading terms largely cancel among each 
other. 

The uncertainty from uncalculated higher order QCD terms is far below 
the foreseeable experimental precision, and the parametric uncertainty from as 
and rrit dominates. An important lesson to be learned from the Zbb vertex in 
one- and two-loop approximation concerns the interplay between "dominant" 
and "sub-leading" pieces: whenever leading and sub-leading terms are of com- 
parable magnitude, inclusion of a sizeable number of terms in the expansion is 
required. The estimate of the final result or of theoretical uncertainties based 
on the first two terms of the series may lead to a wrong or misleading result. 



6 Summary 

Expansion techniques for Feynman amplitudes combined with sophisticated 
computer algebra programs have lead to remarkable progress in multiloop cal- 
culations during the past years. Problems with different mass and energy scales 
can be treated with nested series. Powerful computers allow to evaluate many 
terms in these expansions, and smallness of the expansion parameter is thus 
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no longer required. These techniques have been successfully applied to purely 
hadronic as well as to electroweak observables. 
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